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Liquid  water  management  plays  a  significant  role  in  proton  exchange  membrane  fuel  cell  (PEMFC)  per¬ 
formance,  especially  when  the  PEMFC  is  operating  with  high  current  density.  Therefore,  understanding 
of  liquid  water  behavior  and  flooding  process  is  a  critical  challenge  that  must  be  addressed.  To  overcome 
PEMFC  durability  problems,  a  liquid  water  flooding  process  is  studied  in  the  cathode  side  of  a  PEMFC 
with  straight  parallel  channels  and  a  porous  layer  using  FLUENT®  v6.3.26  software  with  a  volume-of-fluid 
(VOF)  algorithm  and  user-defined-function  (UDF).  The  general  process  of  liquid  water  flooding  within 
this  type  of  PEMFC  cathode  is  investigated  by  analyzing  the  behavior  of  liquid  water  in  porous  layer  and 
gas  flow  channels.  Two  important  phenomena,  the  “first  channel  phenomenon”  and  the  “last  channel 
phenomenon”,  and  their  effects  on  the  flow  distribution  along  different  parallel  channels  are  discussed. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  fuel  cell  is  an  environment  friendly  energy  conversion  device 
that  directly  converts  chemical  energy  into  electrical  energy.  It  is 
regarded  as  a  major  alternative  to  traditional  energy  sources  such 
as  internal  combustion  engines.  Due  to  the  advantages  of  its  low 
operating  temperature  and  high  power  density,  proton  exchange 
membrane  fuel  cell  (PEMFC)  is  considered  one  of  the  potential  can¬ 
didates  for  commercialization. 

For  the  last  two  decades,  although  many  different  types  of 
PEMFC  flow-field  designs  have  been  developed  by  industrial  and 
academic  researchers  [1],  most  of  them  can  be  classified  into  par¬ 
allel,  serpentine,  or  interdigitated  category  [2].  These  different 
designs  were  proposed  to  resolve  the  major  concerns,  such  as  pres¬ 
sure  drop,  uniformity  of  reactant  gas  distribution,  and  liquid  water 
management. 

As  a  critical  research  topic  for  optimizing  the  performance  of 
PEMFCs,  liquid  water  management  is  attracting  great  efforts  from 
many  scientists  and  engineers  [3]. 

Djilali  et  al.  [4]  developed  a  two-dimensional  model  to 
investigate  the  transport  process  in  PEMFCs  to  improve  water  man¬ 
agement.  Previously  reported  analyses  by  Djilali  and  Lu  [5]  also 
studied  the  influence  of  heat  and  water  transport  in  PEMFCs.  A 
three-dimensional  numerical  model  of  straight  gas-flow  channels 
was  studied  by  Dutta  et  al.  [6]  using  the  commercial  computational 
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fluid  dynamics  (CFD)  software  FLUENT®.  Cha  et  al.  [7]  discussed 
the  influence  of  oxygen  concentration  along  gas  flow  channels  on 
PEMFC  performance  by  examining  the  steady-state  gas-flow  phe¬ 
nomena  in  micro  parallel  flow  channels.  Kulikovsky  [8]  numerically 
acquired  the  gas  concentration  of  a  steady-state  flow  along  chan¬ 
nels  in  a  similar  way.  These  studies,  however,  did  not  investigate 
the  influence  of  liquid  water.  Yi  et  al.  [9]  indicated  that  water  vapor 
condensation  is  a  common  phenomenon  on  both  sides  of  a  PEMFC. 
You  and  Liu  [10]  pointed  out  that  a  multi-phase  model  must  be 
employed  in  order  to  obtain  more  practical  results.  Both  of  these 
previous  studies  considered  water  in  vapor  form  [9,10].  Wang  et 
al.  conducted  a  two-phase,  multi-component  mixture  model  of  the 
cathode  of  a  PEMFC  to  deal  with  the  problems  in  a  single-  and  in  a 
two-phase  co-existence  region  [11]  and  introduced  a  technique  to 
measure  in  situ  water  distribution  [12].  Recently,  Dawes  et  al.  [13] 
investigated  the  effects  of  water  with  a  three-dimensional  mixture 
model  of  a  PEMFC,  and  Wang  et  al.  also  developed  a  topologically 
equivalent  pore  network  model  to  study  the  transport  of  liquid 
water  in  PEMFC  [14]. 

Zhou  et  al.  made  significant  efforts  in  investigating  liquid  water 
behavior  in  PEMFCs,  including  in  single  cells  and  in  stacks.  Volume- 
of-fluid  (VOF)  model  was  applied  to  simulate  the  process  of  liquid 
water  removal  using  FLUENT®  [15].  Zhou  et  al.  also  conducted  a 
numerical  simulation  of  two-phase  flow  in  several  different  types  of 
cathode  design,  e.g.,  in  stack  with  parallel  serpentine  channels  [16] 
and  in  stack  with  straight  micro  parallel  channels  [17].  Different 
innovative  gas  diffusion  layers  (GDL),  their  water  removal  charac¬ 
teristics  [18],  the  effects  of  electrode  wettability  on  liquid  water 
behavior  [19],  and  accelerated  numerical  test  of  liquid  behavior 
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Nomenclature 

s  momentum  source  term 

g  acceleration  of  gravity  (ms-2) 

u  velocity  vector  (ms-1) 

I<  surface  curvature 

s  phase  volume  fraction 

p  density 

s  porosity 

/x  dynamic  viscosity  (Pas) 

X  surface  tension  coefficient 

a  gaseous  permeability  (m2) 

Subscripts 
g  gas 

l  liquid 

c  continuity  equation 

m  momentum  equation 


across  GDL  [20]  were  investigated  using  the  VOF  method.  Dji- 
lali  et  al.  [21]  also  used  the  VOF  method  to  study  the  dynamic 
behavior  of  liquid  water  emerging  from  a  GDL  pore  into  a  PEMFC 
gas-flow  channel.  Several  reports  have  contributed  to  the  funda¬ 
mental  understanding  of  liquid  water  behavior  inside  flow  channels 
[16-19],  but  conventional  porous  layers  were  not  included  in  these 
analyses.  Electro-chemical  reactions  have  also  been  neglected. 

A  general  model  for  PEMFCs  [22]  has  been  proposed  and  applied 
to  investigate  two-phase  flow  coupled  with  electro-chemical  reac¬ 
tions,  water  transport  through  the  membrane  and  heat  and  mass 
transfer  for  single  cells  with  serpentine  channels  [22],  single  cells 
with  serpentine-parallel  channels  [23],  and  3-cell  stack  with  par¬ 
allel  design  [24].  These  previous  studies  [22-24]  are  the  state  of 
the  art  for  simulating  PEMFCs;  however,  this  technique  requires 
excessive  computational  time. 

Recently,  Zhou  et  al.  [25]  conducted  numerical  simulation  of 
liquid  water  behavior  in  cathode  of  a  PEMFC  with  serpentine 
channels  and  porous  layer,  using  a  simplified  model  in  which 
electro-chemical  reactions  and  heat  transfer  effects  were  neglected 
to  minimize  the  computation  time.  Also  detailed  experimental  val¬ 
idation  of  this  simplified  cathode  model  was  performed  by  direct 
optical  visualization  method  to  capture  the  motion  and  deforma¬ 
tion  of  liquid  water  with  high  spatial  and  temporal  resolutions.  For 


Fig.  1.  Comparison  of  liquid  water  behaviors  between  numerical  simulation  and  the  experiment  at  t  =  0.400  s  (the  left  panel  shows  the  experimental  setup,  and  the  right 
panel  shows  the  numerical  simulation). 
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Fig.  2.  Simulation  domain. 
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(a)  t  =  0.050  s 

VOF  of  Squid  water 


VOF  of  liquid  water 


VOF  of  liquid  water 


(c)  t  =  0.330s 


(d)  t  =  0.450  s 


VOF  ofliquid  water 


(e)  t  =  0.545  s 


VOF  of  liquid  water 


(f)  t  =  0.590  s 


Fig.  3.  (a-f)  General  process  ofliquid  water  removal:  (a)  t= 0.050  s;  (b)  t  =  0.264  s;  (c)  t= 0.330  s;  (d)  t  =  0.450  s;  (e)  t=  0.545  s;  (f)  t  =  0.590  s. 


reader’s  convenience,  a  sample  comparison  from  this  validation 
is  shown  in  Fig.  1,  which  demonstrates  that  the  results  from  this 
simplified  model  are  in  good  agreement  with  the  corresponding 
experiments. 

From  above  literature  review,  some  key  points  could  be  sum¬ 
marized  as  follows: 


(1 )  Although  significant  efforts  have  been  made  to  understand  the 
liquid  water  effects  on  PEMFC,  the  liquid  water  flooding  process 
still  remain  to  be  investigated  for  all  kinds  of  flow  field  designs. 
For  example,  PEMFC  with  parallel  channel  design  suffers  from 
liquid  water  flooding,  while  there  is  no  literature  available  to 
address  this  important  issue  with  VOF  method. 

(2)  The  VOF  method  with  serpentine  PEMFC  application  has  been 
validated  by  experiments  [25]  and  it  can  be  employed  in  the 
investigation  of  PEMFC  with  other  different  types  of  flow  field 
designs. 

(3)  There  are  three  different  ways  to  introduce  liquid  water  in  VOF 
modeling  of  PEMFC:  presetting  certain  amounts  of  droplets 


[15-17,22,23,25]  or  certain  liquid  water  film(s)  [16,17],  inject¬ 
ing  liquid  water  droplets  through  a  pinhole  [21,24],  and  setting 
liquid  water  flowrate  using  the  accelerated  numerical  method 
[20].  The  third  one  can  provide  significant  insights  into  the  fun¬ 
damental  understanding  ofliquid  water  flooding  process  inside 
PEMFC  with  an  acceptable  computation  time. 

Using  the  numerical  methodology  developed  and  validated 
in  the  present  authors’  previous  works  [20,25]  as  an  engineer¬ 
ing  research  tool,  the  present  paper  investigates  the  liquid  water 
behavior  inside  a  PEMFC  cathode  with  straight  parallel  chan¬ 
nels  and  a  porous  layer  to  gain  the  fundamental  understanding 
on  the  liquid  water  flooding  process  within  this  type  of  PEMFC 
design. 

In  the  following  sections,  the  numerical  simulation  domain, 
mathematical  model,  mesh  independency,  and  initial  settings  are 
introduced.  These  sections  are  followed  by  a  detailed  discussion  of 
liquid  water  flooding  process,  porosity  effects  on  drainage  perfor¬ 
mance,  and  conclusions. 


X.  Wang,  B.  Zhou /Journal  of  Power  Sources  196(2011)  1776-1794 


1779 


VOF  of  liquid  water 


X,  m 


VOF  of  liquid  water 


0.95 

0.85 

0.75 

0.65 

0.55 

0.45 

0.35 

0.25 

0.15 

0.05 


(a)  t  =  0.200  s 
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Fig.  4.  (a-d)  Liquid  water  distribution  at  Z= 0.00029  m:  (a)  t=  0.200  s;  (b)  t=  0.251  s;  (c)  t=  0.284  s;  (d)  t=  0.304  s. 


2.  Numerical  model  setup 

2A.  Numerical  simulation  domain  and  boundary  conditions 

The  numerical  simulation  domain  represents  as  full-scale  cath¬ 
ode  side  geometry  of  a  single  PEMFC.  Fig.  2  illustrates  the  schematic 
of  the  numerical  simulation  domain,  which  contains  gas  flow  chan¬ 
nels  and  a  porous  layer.  The  channels  are  in  a  parallel  pattern,  which 
permits  a  low  overall  pressure  drop  between  the  gas  inlet  and  the 
outlet.  There  are  seven  parallel  channels,  with  a  0.001  m  distance 


between  adjacent  channels.  The  length  of  the  cathode  from  the  air 
inlet  to  the  outlet  is  0.06  m  in  the  X-direction.  The  width  of  the 
gas  flow  channel  domain  is  0.02  m  in  the  Y-direction.  The  cross- 
section  of  the  channels  is  0.002  m  (in  the  X-direction)  x  0.002  m 
(in  the  Y-direction)  x  0.0017  m  (in  the  Z-direction).  The  dimen¬ 
sions  of  the  porous  layer  are  0.024  m  (in  the  X-direction)  x  0.024  m 
(in  the  Y-direction)  x  0.0003  m  (in  the  Z-direction);  these  are 
the  typical  dimensions  of  a  fuel  cell  component.  A  liquid  water 
injection  channel  (LWIC)  with  dimensions  of  0.0004  m  (in  the 
X-direction)  x  0.0004  m  (in  the  Y-direction)  x  0.000725  m  (in  the 


0  006 


X,  m 


8 


0024 


VOF  of  liquid  water 


- 

I  0.95 
085 

075 

065 

055 

045 

035 

025 

015 

005 

Pressure 


^  210 
—  170 
—\  130 

■  90 
50 


(a)  t  =  0.200  s 


(b)  t  =  0.200  s 


Fig.  5.  (a)  VOF  and  velocity  distribution  at  Z= 0.00029  m  and  (b)  pressure  distribution  at  Z= 0.00029  m.  (a)  t  =  0.200  s  and  (b)  t  =  0.200  s. 
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Z-direction)  is  positioned  at  the  inlet  joint  (defined  in  Fig.  2).  This 
LWIC  was  used  to  simulate  liquid  water  injection  effects  in  previ¬ 
ous  studies  [23]  from  Zhou’s  group.  In  the  present  study,  the  LWIC 
is  blocked  because  there  no  liquid  water  is  injected  through  it  in 
this  study,  and  this  geometry  is  kept  so  that  the  results  for  the 
current  study  (where  the  liquid  water  is  supplied  from  the  back 
surface  of  the  porous  layer  as  shown  in  Fig.  2)  can  be  compared 
to  the  results  of  previous  studies  (where  liquid-water  injection  is 
performed  through  the  LWIC).  The  direction  of  gravity  is  along  the 
negative  Z-direction,  which  is  the  normal  direction  of  the  porous 
layer,  as  shown  in  Fig.  2. 

A  no-slip  boundary  condition  is  applied  to  the  walls.  The  air 
inlet  is  set  with  a  flowrate  of  2.00  x  10-5  kgs-1  (the  flow  is  nor¬ 
mal  to  the  inlet  boundary)  and  the  liquid  water  inlet  (back  surface 
of  the  porous  layer  shown  in  Fig.  2)  is  set  with  a  flowrate  of 
1.70  x  10-4  kgs-1.  The  pressure  outlet  is  applied  at  the  outlet 
boundary.  The  contact  angles,  which  are  set  in  the  side  walls  of 
channels,  upper  walls  of  channels  and  the  rest  of  walls,  are  40°,  43° 


and  140°,  respectively.  The  porosities  of  the  porous  layers  used  in 
this  model  are  0.3,  0.5  and  0.8  for  the  three  different  cases  con¬ 
ducted  in  the  present  study. 

2.2.  Computational  methodology 

The  simulation  of  a  three-dimensional,  unsteady,  two-phase 
flow,  in  which  the  air  is  treated  as  gaseous  phase  and  the  liquid 
water  is  treated  as  liquid  phase  (the  two  phases  are  assumed  to 
be  immiscible  in  this  model),  is  conducted  using  FLUENT®  v6.3.26 
software.  The  VOF  method  is  applied  to  track  the  air-water  flow 
interface  inside  the  computational  domain,  and  the  UDF  is  applied 
to  introduce  the  source  term. 

In  this  model,  the  fluid  flow  is  laminar  because  of  the  low 
flow  velocity  and  the  small  size  of  the  channels.  In  respect 
that  heat  transfer  and  phase  change  are  not  considered  in  this 
study  in  order  to  save  the  computation  time,  the  unsteady,  lam¬ 
inar  flow  model  is  governed  by  the  following  continuity  (mass) 
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Fig.  6.  (a-h)  General  process  by  which  liquid  water  emerges  into  the  channels:  (a)  t=  0.200  s;  (b)  t  =  0.251  s;  (c)  t  =  0.284s;  (d)  t= 0.304  s;  (e)  t=  0.322  s;  (f)  t= 0.388  s;  (g) 
t  =  0.500s;  (h)t  =  0.655s. 
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Fig.  6.  ( Continued ) 


and  momentum  equations  to  describe  the  fluid  flow  transport 
process. 

The  continuity  (mass)  equation  is  expressed  as 
^+V.(epii)  =  Sc  (1) 


where  s  is  the  porosity  of  the  porous  media  (in  the  gas  flow  chan¬ 
nels  s  =  1 ).  The  first  term  on  the  left  hand  side  is  the  transient  term 
standing  for  the  change  of  mass  with  time,  the  second  term  is  the 
convection  term  standing  for  the  change  in  mass  flux,  and  the  term 
on  the  right  hand  side  is  the  mass  source  term  which  equals  to  zero 
because  phase  changes  are  not  taken  into  account  in  this  model. 


In  the  VOF  model,  the  volume  fractions  of  gas  and  liquid  water 
are  defined  as  s/  and  sg,  respectively,  and 

si  +  sg  —  1  (2) 

The  variables  and  properties  are  shared  by  both  phases  and  are 
calculated  by  a  volume-averaged  method.  The  average  density  and 
viscosity  for  each  computational  cell  which  are  between  those  of 
gaseous  and  liquid  phases,  are  defined  as 

p  =  s,p,  +  sgpg  (3) 

and 

M  =  $i!M  +  ssp-g  (4) 
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Fig.  6.  ( Continued ) 


In  this  VOF  model,  the  interface  between  gaseous  and  liquid 
phases  is  tracked  by  solving  a  continuity  equation,  i.e.,  Eq.  (5),  for 
the  volume  fraction  of  liquid  water,  which  is  solved  over  the  entire 
computational  domain.  Then  the  volume  fraction  of  gas  is  com¬ 
puted  based  on  Eq.  (2). 


The  Geo-Reconstruct  scheme  is  applied  to  smooth  the  interface 
between  two  phases  using  a  piecewise-linear  approach,  which  is 
more  accurate  than  other  methods  in  FLUENT®  [26]. 

A  single  momentum  equation,  i.e.,  Eq.  (6),  is  implemented  in  the 
VOF  model,  and  the  gaseous  phase  and  the  liquid  phase  share  the 
same  resulting  velocity  field. 


d(esiPi) 


+  V  ■  ( ssipiu )  =  0 


(5) 


9 

— (^pu)  +  V  •  ( spuu )  =  -£Vp  +  V[s/lVu]  +  £Sm 


(6) 


dt 


X.  Wang,  B.  Zhou/Journal  of  Power  Sources  196(2011)  1776-1794 


1783 


0  0.006  0.012  0.018  '  0.024 

X,  m 


1  m/s  Y  -  0.005  m 


parallel  channel  #4 
X  :  Z  -  1  : 1 


0.2  m/s 
0.002  m/s 


X,  m 


(d)  t  =  0.304  s 

Fig.  6.  ( Continued ) 


In  this  momentum  equation,  the  terms  on  the  left  hand  side 
represents  the  change  of  momentum  with  time  and  the  advection 
momentum  flux,  respectively,  while  the  terms  on  the  right  hand 
side  are  the  momentum  due  to  pressure,  viscosity  and  sources. 

In  the  flow  channels  of  this  model,  the  surface  tension  force,  as 
well  as  gravity  force,  is  coupled  to  the  momentum  equation  by  the 
source  term  Sm  using  the  continuum  surface  force  (CSF)  method 
[27]  which  helps  smooth  the  interface  in  the  VOF  model,  expressed 
as 


In  the  porous  layer,  one  more  term  -(/z/a)u,  where  a  stands  for 
the  permeability  of  the  porous  layer,  is  added  in  the  source  term 
because  of  the  viscous  effect  though  the  porous  layer,  or  the  Darcy 
force.  So  the  source  term  for  the  momentum  equation  in  the  porous 
layer  is 


Sm 


—  Pg  +  XK 


pVsi 

Pi  +  Pg/2 


(8) 


2.3.  Grid  independency 


Sm 


—  Pg  +  XK 


pVsi 

Pi  +  Pg/2 


(7) 


where  x  is  the  surface  tension  coefficient  and  I<  is  the  curvature. 


The  computation  domain  is  meshed  into  290,352  cells  with  a 
maximum  volume  of  5.54861  lxl  0-12  m3  and  a  minimum  volume 
of  5.0  x  10-13m3  for  each  cell.  The  grid  size  measures  approxi¬ 
mately  8.3  x  10-4  m  in  both  the  X-  and  Y-directions,  whereas  the 
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dimensions  in  the  Z-direction  are  1.4x1 0-4  m  in  the  gas  flow  chan¬ 
nel  domain  and  5  x  10-5  m  in  the  porous  layer  domain.  Zhou  et  al. 
[17]  implemented  a  grid  check  method  by  increasing  and  decreas¬ 
ing  the  number  of  grid  cells  by  certain  percentages.  And  this  grid 
size  has  been  validated  for  PEMFCs  by  the  comparison  of  numer¬ 
ical  results  with  experimental  results  from  the  author’s  previous 
studies  [25]. 


3.  Results  and  discussion 

In  this  study,  three  cases  with  different  porosities  (0.3,  0.5  and 
0.8)  have  been  investigated.  The  general  water  flooding  processes 
between  these  three  cases  are  quite  similar;  therefore,  in  the  fol¬ 
lowing  sections,  the  model  with  a  porosity  of  0.3  is  used  for  the 
general  discussion.  The  porosity  effect  on  flooding  process  is  also 
addressed. 
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3. 1.  General  process  of  liquid  water  removal 

Fig.  3(a)-(f)  shows  the  main  process  of  water  movement  over 
time  for  some  selected  time  instants.  The  general  process  of  liquid 
water  removal  can  be  divided  into  the  following  sub-processes: 


(1)  Liquid  water  is  supplied  from  the  back  surface  of  the  porous 
layer  by  a  constant  flowrate  to  simulate  the  water  generation 
through  electro-chemical  reaction  (Fig.  3(a)). 


(2)  This  supplied  liquid  water  first  enters  the  porous  layer  and 
moves  around  with  in  the  porous  layer  (Fig.  3(b)). 

(3)  Liquid  water  from  the  porous  layer  emerges  into  gas  flow  chan¬ 
nels  from  the  peripheral  edges  of  lands  (solid  between  adjacent 
parallel  channels  defined  in  Fig.  2)  and  the  peripheral  edges  of 
the  frames  (Fig.  3(c)). 

(4)  The  liquid  water  emerging  into  the  channels  accumulates  at  the 
outlet  joint  (Fig.  3(d)). 

(5)  The  liquid  water  at  outlet  joint  starts  to  move  into  the  outlet 
channel  (Fig.  3(e)). 
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(6)  Liquid  water  drains  out  through  the  outlet  channel 
(Fig.  3(f)). 

3.2.  Liquid  water  emerging  process 

3.2  A.  Liquid  water  emerging  from  the  peripheral  areas  of  lands 
and  frames 

As  shown  in  Figs.  4(a)-(d)  and  5(a)  and  (b),  a  plane  from 
the  porous  layer  domain  (with  projective  channel  edges)  at 
Z= 0.00029  m,  which  is  very  close  to  the  interface  between  porous 
layer  and  channels,  is  extracted  to  explain  the  main  process  of 
liquid  water  entering  into  the  channels  from  the  porous  layer. 
Fig.  4(a)  shows  that  the  water  initially  emerges  from  the  outskirts 


of  the  porous  layer  (the  porous  layer  area  under  the  frames)  and 
from  the  left-bottom  and  right-top  corners  of  the  areas  near  the 
flow  channels.  Fig.  5(a)  and  (b)  gives  the  velocity  and  pressure 
distributions  at  the  same  time  instant  that  is  used  in  Fig.  4(a) 
(t= 0.200  s).  Fig.  5(a)  shows  that  liquid  water  accumulates  in  the 
outskirts  and  corners  of  the  porous  layer  due  to  the  constraint  of 
the  porous  layer  domain,  so  the  pressure  in  these  areas  is  higher 
than  that  in  other  areas  of  the  porous  layer  (Fig.  5(b)),  resulting  in 
a  flow  from  the  outskirt  areas  to  the  peripheral  edges  of  the  frames 
(Fig.  5(a))  inside  the  porous  layer.  This  flow  explains  why  liquid 
water  emerges  initially  from  the  peripheral  edges  of  the  frames,  as 
shown  in  Fig.  4(b).  Fig.  4(b)-(d)  shows  that  water  enters  into  the 
channels  from  the  edges  of  lands  (the  ribs  between  two  adjacent 
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flow  channels  as  defined  in  Fig.  2)  and  fills  in  the  areas  underneath 
the  channels. 

The  liquid  water  emergence  from  the  porous  layer  into  chan¬ 
nels  can  be  explained  by  Fig.  6(a)-(h),  which  shows  the  velocity 
vectors  and  VOF  on  the  X-Z  planes  located  at  Y=  0.01 9  m,  0.01 2  m, 
and  0.005  m  with  different  reference  velocity  vectors. 

The  liquid  water  from  the  porous  layer  enters  into  the  channels 
from  the  edges  of  the  lands  because  the  liquid  water  generation  rate 
is  uniform  on  the  back  surface  of  the  porous  layer,  but  the  liquid 
water  underneath  the  lands  cannot  directly  enter  into  the  channels 
in  the  manner  that  the  liquid  water  underneath  the  middle  of  the 
channels  enters.  Actually,  the  liquid  water  underneath  the  lands 
has  to  detour  along  the  interface  between  the  lands  and  the  porous 


layer  then  into  the  channels  (Fig.  6(a)-(d))  because  the  lands  are 
solid  and  only  allow  the  conduction  of  electrons.  The  liquid  water 
emerges  into  the  channel  directly  through  the  porous  layer  under¬ 
neath  the  middle  of  the  channels  (Fig.  6(e)),  moving  with  the  main 
gas  flow  along  the  parallel  channels,  and  gradually  merges  with 
the  liquid  water  that  has  accumulated  around  the  corners  of  the 
channels  (Fig.  6(f)-(h))  due  to  surface  tension  and  wall  adhesion 
effects. 

The  emergence  of  liquid  water  from  the  porous  layer  into  the 
channels  has  direct  effects  on  the  gas  flow  inside  the  channels,  in 
the  primary  flow  direction  (the  main  gas  flow  along  the  channels) 
by  blocking  the  channels,  and  on  the  cross-sections  by  generating 
various  secondary  flows,  as  shown  in  Fig.  6(c)-(h). 
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3.2.2.  Dean  vortex  evolution 

For  the  X-Z  plane  at  Y=  0.01 9  m,  i.e.,  the  plane  cut  near  the 
entrance-turning-area  shown  in  Fig.  6(a)-(h),  it  can  be  noticed  that 
there  are  Dean  vortices  in  the  parallel  channels.  The  lower  vortex  of 
the  vortices  gradually  shrinks  as  the  liquid  water  emerges  from  the 
porous  layer  into  the  channels.  From  the  first  to  the  last  channels, 
the  strength  of  the  vortices  increases  because  the  maldistribution 
of  the  gas  flow  among  the  parallel  channels  increases;  this  phe¬ 
nomenon  will  be  discussed  in  the  later  section  titled  “Mechanism 
of  first  and  last  channel  phenomena”. 

For  the  X-Z  plane  at  Y=  0.01 2  m,  i.e.,  the  plane  cut  in  the  middle 
of  the  domain  in  the  Y-direction,  there  are  distorted  Dean  vortices 
in  the  last  channel.  The  lower  vortex  also  gradually  shrinks  as  the 
liquid  water  emerges  from  the  porous  layer  into  the  channel.  At 
a  certain  time  (Fig.  6(h)),  the  lower  vortex  disappears  completely, 


which  shows  the  effects  of  gas-liquid  interactions.  There  are  no 
(distorted)  Dean  vortices  from  parallel  channel  #1  to  #6  because 
these  areas  are  removed  from  the  influence  of  the  Dean  vortices 
around  the  entrance-turning-areas. 

For  the  X-Z  plane  at  Y=  0.005  m,  i.e.,  the  plane  cut  near  the  exit¬ 
turning-area,  there  are  no  Dean  vortices  because  this  plane  is  too 
far  away  from  the  entrance-turning-area. 

3.3.  General  liquid  water  removal  process  in  gas  flow  channels 

The  liquid  water  removal  process  in  gas  flow  channels  can  be 
explained  by  Fig.  7.  The  X-Y  view  of  the  VOF  of  the  liquid  water 
inside  the  channels  is  shown  in  Fig.  7(a-l)-(n-l),  and  the  water 
distributions  in  the  three  selected  cross-sections  with  Y=  0.005  m, 
Y=  0.01 2  m  and  Y  =  0.01 9  m  are  defined  in  Fig.  7(a-l )  and  shown  in 
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Fig.  7.  (a-n)  General  process  of  liquid  water  removal  in  channels:  (a-1)  t= 0.500  s,  (a-2)  t  =  0.500  s;  (b-1)  t=  0.523  s,  (b-2)  t  =  0.523  s;  (c-1)  t= 0.544  s,  (c-2)  t  =  0.544  s;  (d-1) 
t  =  0.546s,  (d-2)  t  =  0.546s;  (e-1)  t  =  0.555s,  (e-2)  t=0.555s;  (f-1)  t  =  0.590s,  (f-2)  t=0.590s;  (g-1)  t  =  0.655s,  (g-2)  t  =  0.655s;  (h-1)  t=0.724s,  (h-2)  t  =  0.724s;  (i-1)  t  =  0.759s, 
(i-2)  t  =  0.759  s;  (j-1)  t= 0.920  s,  (j-2)  t  =  0.920s;  (k-1)  t  =  0.965  s,  (k-2)  t  =  0.965  s;  (1-1)  t  =  0.994s,  (1-2)  t= 0.994  s;  (m-1)  t=  1.171s,  (m-2)  t=  1.171s;  (n-1)  t=  1.300  s,  (n-2) 
t=  1.300  s. 
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the  right  column  of  Fig.  7(a-2)-(n-2);  this  information  will  help  the 
reader  understand  the  status  of  the  liquid  water  blockage  across 
parallel  channels. 

The  liquid  water  coming  from  the  porous  layer  near  the  periph¬ 
eral  zones  of  the  four  frames  and  the  lands  shown  in  Fig.  7(a) 
gradually  develops  into  two  liquid  water  streams,  with  the  first 
stream  accumulating  at  the  bottom  of  the  outlet  manifold  and 
the  second  stream  accumulating  on  the  side  wall  of  the  last  par¬ 
allel  channel  (#7).  These  two  streams  meet  near  the  outlet  joint 


(Fig.  7(a))  and  gradually  form  a  wheel  shaft  shape  (two  wheels 
connected  by  one  shaft  in  the  middle)  (Fig.  7(b)).  This  shape  is 
formed  due  to  the  combined  effects  of  wall  adhesion  and  liquid 
water  surface  tension.  At  the  same  time,  the  air  streams  from  all 
directions  also  meet  and  merge  into  one  stream  at  the  outlet  joint; 
this  combined  air  stream  continuously  strikes  the  left  side  of  the 
wheel  shaft,  forcing  the  wheel  shaft  to  move  along  the  outlet  chan¬ 
nel  and  break  into  a  pair  of  wheel  shafts  (Fig.  7(c)  and  (d)).  Then,  this 
wheel  shaft  shape  deforms  into  a  top-bottom  stratified  flow  pat- 


1790 


X.  Wang,  B.  Zhou  /  Journal  of  Power  Sources  196  (201 1)  1 776-1794 


(h-1)t  =  0.724  s 


VOF  of  liquid  water 


VOF  of  liquid  water 


VOF  of  liquid  water 


(i-1)t  =  0.759  s 


0-1)  t  =  0.920  s 


0-2)  t  =  0.920  s 


Fig.  7.  ( Continued ) 


tern  while  moving  along  the  outlet  channel  to  drain  (Fig.  7(e)  and 
(f)).  These  processes  repeat  for  a  period  of  time  (Fig.  7(g)-(m)),  and 
the  liquid  water  eventually  moves  along  the  outlet  channel  with  a 
top-stratified  pattern  (Fig.  7(n)).  Within  this  time  period,  two  inter¬ 
esting  phenomena  occur;  i.e.,  the  first  and  the  last  parallel  channel 
phenomena.  These  phenomena  will  be  explained  in  the  following 
section. 

3.3.1.  First  and  last  parallel  channel  phenomena 

At  a  certain  time  instant  (t  =  0.590  s  for  the  present  case)  and 
near  the  exit-turning-area  (defined  in  Fig.  2)  of  the  last  parallel 
channel,  the  liquid  water  streams  coming  from  the  peripheral  edge 
of  the  land  and  the  edge  of  the  right  frame  merge  together  (Fig.  7(f)). 
Then,  as  shown  in  Fig.  7(g),  a  similar  merging  process  happens  near 
the  entrance-turning-area  (defined  in  Fig.  2)  of  the  first  parallel 
channel  due  to  the  vortex  that  is  caused  by  the  inertia  of  the  inlet 
air  flow  (Fig.  8(a)). 

Almost  immediately  after  the  blockage  of  the  first  parallel  chan¬ 
nel  occurs  near  the  entrance-turning-area,  liquid  water  starts  to 


accumulate  around  the  left-bottom  corner  of  the  outlet  manifold.  At 
the  moment  when  this  water  front  touches  the  left-bottom  corner 
of  the  first  land,  the  first  parallel  channel  is  blocked  from  the  bot¬ 
tom  (Fig.  8(b)).  Then,  this  liquid  water  front  continuously  expands 
along  the  first  channel  from  the  bottom  up  (Fig.  7(k)-(n))  due  to  the 
lower  pressure  zone  caused  by  the  blockage  near  the  first  entrance¬ 
turning-area  that  is  established  just  above  the  water  front  inside  the 
first  channel. 

Around  the  time  instant  t= 0.724  s  (Figs.  7(h)  and  9(a))  at  the 
entrance-turning-area  of  the  last  parallel  channel,  the  liquid  water 
stream  coming  from  the  peripheral  edge  of  the  land  and  the  stream 
emerging  from  the  edge  of  the  right  frame  merge  together.  Almost 
immediately  after  this  merge,  the  accumulated  water  in  the  right- 
top  corner  of  channels  starts  to  splash  into  small  droplets  as  it  is 
struck  by  the  gas  flow  from  the  inlet  manifold  (Fig.  7(h)-(k)).  There 
is  a  low  pressure  zone  around  the  entrance-turning-area  of  the 
last  parallel  channel,  which  attracts  droplets  of  liquid  water  from 
around  the  corner  into  the  main  flow  stream.  These  small  droplets 
continue  to  move  toward  the  land  wall  and  adhere  to  it  due  to 
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the  wall  adhesion  effect.  The  interaction  between  the  airflow  and 
the  liquid  water  in  the  vicinity  of  corners  can  provide  strong  shear 
stress  and  can  fragment  some  of  the  liquid  water  on  the  gas-liquid 
water  interface  into  small  droplets.  These  small  droplets  flow  along 
the  side  wall  of  the  right  frame  due  to  inertial  forces  and  wall  adhe¬ 
sion,  and  then,  the  small  liquid  water  droplets  merge  into  large 
droplets  (Figs.  7(1)  and  9(b))  that  continuously  flow  along  the  side 
wall  of  the  right  frame,  reaching  the  corner  near  the  outlet  joint 
and  accumulating  there.  This  accumulation  is  because  that  near  the 
joint  the  surface  tension  forces  and  wall  adhesion  forces  are  large 
due  to  the  large  curvature  of  the  wall. 

Around  time  t=  1.171  s,  as  shown  in  Fig.  7(m),  the  liquid  water 
that  is  accumulating  around  the  corner  of  the  right  frame  starts  to 
drain  out  along  the  top  side  wall  of  the  outlet  channel.  This  type 
of  flow  is  known  as  top-stratified  drainage  flow,  and  this  effect 
becomes  dominant  after  this  point  in  time. 

3.3.2.  Mechanism  of  first  and  last  channel  phenomena 

Fig.  9  shows  the  VOF  and  velocity  distribution  on  three  planes 
( Y  =  0.005  m,  Y=  0.01 2  m  and  Y=  0.01 9  m).  At  t = 0.200  s  (Fig.  1 0(a)), 
when  all  liquid  water  still  remains  in  the  porous  layer,  the  air  distri¬ 
bution  among  parallel  channels  increases  with  increasing  distance 
from  the  inlet  joint  (defined  in  Fig.  2),  with  the  last  channel  having 
the  maximum  flowrate  and  velocity  and  the  first  channel  having 
the  minimum  flowrate  and  velocity.  Fundamentally,  this  maldis¬ 


tribution  of  gas  flow  is  caused  by  the  parallel  channel  design;  this 
point  is  also  discovered  by  Wang  et  al.  in  his  recent  publication 
[28].  The  maldistribution  of  reactant  gas  (air)  flow  then  creates  a 
maldistribution  of  liquid  water  in  the  parallel  channels  and  porous 
layers.  The  maldistribution  of  liquid  water,  in  return,  further  affects 
the  of  the  reactant  gas  flow  (Fig.  10(b)).  This  emphasizes  that  liquid 
water  management  is  still  one  of  the  critical  challenges  for  PEMFC 
designs. 

3.4.  Porosity  effect  on  drainage  performance 

To  investigate  the  influence  of  the  porosity  of  the  porous  layer 
on  the  water  flooding  process,  three  cases  with  different  porosities 
(0.3,  0.5  and  0.8)  were  studied.  The  general  liquid  water  flood¬ 
ing  processes  are  quite  similar  among  these  three  cases,  whereas 
the  drainage  performance  parameters,  such  as  the  water  emerg¬ 
ing  time  (the  instant  when  the  liquid  water  from  the  porous  layer 
emerges  into  the  gas  flow  channels)  and  the  water  departure  time 
(the  instant  when  the  liquid  water  begins  to  depart  from  the  outlet 
of  the  fuel  cell),  are  different.  The  water  emerging  time,  departure 
time  and  porosity  for  these  three  cases  are  shown  in  Fig.  11.  The 
results  reveal  that  both  times  change  in  and  approximately  linear 
fashion  with  respect  to  porosity,  i.e.,  larger  porosities  give  longer 
times  because  the  porous  layer  with  larger  porosity  can  accommo¬ 
date  more  liquid  water  for  a  longer  time. 
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Fig.  9.  (a,  b)  Enlarged  view  of  liquid  water  splashing  phenomenon  at  right-top  corner  of  channels:  (a)  t=  0.724  s  and  (b)  t  =  0.994  s. 
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Fig.  10.  (a,  b)  VOF  and  velocity  distribution  on  three  planes  (7=0.005  m,  7=0.012  m  and  7=  0.019  m):  (a)  t  =  0.200  s  and  (b)  t=  1.171  s. 
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Fig.  11.  Porosity  affection  on  drainage  performance. 


4.  Conclusions 

The  liquid  water  flooding  process  in  the  porous  layer  and  in  the 
parallel  channels  on  the  cathode  side  of  a  PEMFC  was  investigated 
by  conducting  a  three-dimensional,  unsteady  two-phase  flow  using 
a  VOF  model  in  FLUENT®.  The  porosity  effect  on  the  liquid  water 


draining  process  was  also  studied.  The  main  conclusions  can  be 
summarized  as  follows: 

(1 )  Liquid  water  flooding  behavior  in  the  cathode  is  directly  influ¬ 
enced  by  the  parallel  channel  design. 

(2)  The  liquid  water  accumulation  usually  starts  from  the  zones 
under  the  frames  inside  the  porous  layer. 

(3)  The  liquid  water  from  the  porous  layer  first  enters  into  the  air 
flow  channels  through  the  edges  of  lands  and  frames. 

(4)  There  are  two  important  phenomena  (the  first  channel  phe¬ 
nomenon  and  last  channel  phenomenon)  that  significantly 
affect  the  liquid-water  flow  patterns  in  the  outlet  channel. 

(5)  The  porosity  of  the  porous  layer  is  one  of  the  important  fac¬ 
tors  which  affect  the  water  removal  process.  The  porosity  can 
change  the  drainage  performance  by  changing  the  time  for  liq¬ 
uid  water  to  enter  into  and  drain  out  of  the  channels. 

In  the  future,  the  liquid  water  flooding  processes  for  cathodes 
with  different  channel  designs  will  be  investigated  and  compared 
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to  this  parallel  channel  study.  This  comparison  will  aid  in  gaining  a 
fundamental  understanding  of  the  effects  of  channel  design  on  the 
liquid  water  removal  process. 
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